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'do'- ABSTRACT 

o ■ 

O 

^ ■ We compute the mid-infrared extinction law from 3.6 — 24/im in three molecu- 
lar clouds: Ophiuchus, Perseus, and Serpens, by combining data from the "Cores 

^ . to Disks" Spitzer Legacy Science program with deep JHKg imaging. Using a 

■ new technique, we are able to calculate the line-of-sight extinction law towards 

^ . each background star in our fields. With these line-of-sight measurements, we 

Q^. create, for the first time, maps of the deviation of the data from two extinc- 

Q . tion law models. Because our maps have the same spatial resolution as our 

^ . tion of the total column density. In the Spitzer IRAC bands, 3.6 — 8 /im, we see 

evidence for grain growth. Be low Ak^ = 0.5, our extinction law is well-fit by the 



Weingartner fc Draind (1200 ll ) Ry = 3.1 diffuse interstellar medium dust model. 



^ . As the extinction increases, our law gradually flattens, and for A^^ > 1, the data 

are more consistent with the Weingartner & Draine Ry = 5.5 model that uses 
larger maximum dust grain sizes. At 24 /xm, our extinction law is 2 — 4x higher 
^ ! than the values predicted by theoretical dust mode ls, but is more consistent with 

OO ! the observational results of iFlaherty et al.l (120071 ) . Lastly, from our maps 

^ ' we identify a region in Perseus where the IRAC extinction law is anomalously 

high considering its column density. A steeper near-infrared extinction law than 
the one we have assumed may partially explain the IRAC extinction law in this 
region. 
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Introduction 



Dust is a ubiquitous component of the interstellar medium. It is important to under- 
stand the dust because it attenuates, or extincts, the light from background objects seen 
through the dust. The amount of extinction varies with both wavelength and position. 
Because of the constancy o f the gas-to-dust r atio, the quantity of dus t gives you a direct 
measure of the cloud's mass (IBohlin et al.lll978l ). Note, however that the lBohlin et al.l (119781 ) 
equation does depend on the E{B — V) color excess, which can potentially vary depending 
on interstellar medium properties such as metallicity. The dust extinction must also be ac- 
counted for when computing the luminosities for protostars within molecular clouds. The 
wavelength dependence of the extinction, known as the extinction law, directly relates to 
dust grain properties such as size, composition, and the presence or lack of icy mantles on 
the grains. 



Many authors have found that the near-infrared extinction law has t he form of a p ower 
law. Ax oc A"'^, between A ~ 1 /im and ~ 4 /xm with /5 = 1.6 — 1.8 (jPrainel l2003l . and 



references th erein). The extinc t ion la w in the mid-infrared wavelengths has not been as 
well studied. iRieke fc LebofskyI (119851 ) found the near-infrared power law extended out to 
~ 7 fim before the extinctions increased again towards the 1 /im s ilicate absorption peak. 



More recent results have found a different behavior. iLuta (119991) used ISO t o me asure 
hydrogen recombinat ion lines towards the Galactic center and llndebetouw et al.l (120051 ) and 
Flaherty et al.l (120071 ) used the Spitzer Space telescope to measure the mid-infrared extinction 
law in several regions. These authors all find a s imilar extinction law in the mid-infrared 
which is much flatter than the lRieke fc LebofskyI (119851 ) result. 



To date no one has studied the changes in the extinction law within a large region with 
high spatial resolution. In this paper we use Spitzer and deep ground-based JHKg obser- 
vations to probe changes in the dust properties within three molecular clouds: Ophiuchus, 
Perseus, and Serpens. We compare our results with different dust models to investigate the 
relationship between dust grain properties and total extinction (column density). In §[2] we 
describe our observations. Then, in §[3] we discuss our data reduction method, including how 
we identified background stars and computed line-of-sight (LOS) extinctions. In §11115] we 
create extinction and maps and compare them to each other. We also compute the aver- 
age extinction law from 3.6 — 24 fim for different ranges of extinction. Finally, we summarize 
our results in §[3 
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Observations 



We mapped regions witliin tlie Opliiuclius, Perseus, and Serpens molecular clouds in 
eight wavebands from 1.25 — 24 /im. Our Spitzer observations (3.6 — 24 /im) were part 
of the Spitzer Legacy Sc ience program "From Molecular Cores to Planet- Forming Disks'' 



c2d) jEvans_eial.ll2003h. The Spitzer o bservations have been previously publish ed for 



Serpens ( iHarvey et al. 



2006 



Ophiuch us flPadgett et all 120071). P erseus (IJgirgensen et all l2006l : iRebull et all 120071 ). and 



2007bl ). To complement these data, we took JHKg observa- 



tions at Kitt Peak National Observatory in three observing runs from September 2004- June 
2006. We used th e FLoridA M ulti-object Imaging Near-ir Grism Observations Spectrometer 
(FLAMINGOS) JElstonlllQQs] ) on the 4-meter telescope. 



The c2d observations of Ophiuchus, Perseus, and Serpens cover approximately 11 deg^ 
with both Spitzer^s IRAC (3.6 — 8 fim) and MIPS (24 /xm) instruments. Such a large area 
could not feasibly be observed by FLAMINGOS. Instead, we focused on selected regions 
within these clouds. Our goal was to explore changes in the extinction law within molecular 
clouds. Therefore, we chose contiguous regions with a range of extinction values, from low 
to high. Furthermore, we tried to include some star-forming regions within our mapped 
areas since these regions can alter the dust properties as well. In Table [1] we list the basic 
properties of these clouds including the total area mapped with our JHKg observations. 

In Figures [1] - [3] we plotted a three color image of each region using the c2d data. 
The red, green, and blue emission correspond to the 24, 8, and 3.6 /im channels on Spitzer, 
respectively. The green outline denotes the area we mapped in the JHKg bands. 

In Ophiuchus we mapped a ~ 10' x 60' region to the northeast of L1688. The western 
side of our map, near L1688, contains many young stellar objects (YSOs) and several mid- 
infrared dark cores are seen in the middle and on the eastern side. The 24 /im emission is 
particularly bright on the western half of our map causing a saturated region to appear cyan 
in our color image. 

Our Perseus observations also have a significant amount of dust nebulosity. A partic- 
ularly bright region near the center of our map contains the IRAS source 03382+3145. To 
the southeast of this region is a cluster of bright red sources. Based on these sources' red 
appearance in Figure they likely have an infrared excess suggesting they are young stellar 
objects. These sources were not resolved by IRAS, but IRA S 03388-f3139, a known infrared 
source is located at this position. This cluster is identified by |j0rgensen et al.l (120061 ). Several 
other isolated infrared excess sources exist within the field. Just to the east of our observed 
region is the well-known star-forming region IC348. 



Finally, in Serpens we mapped a region surrounding "Cluster B" , a group of YSOs to 
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the south of the Serpens core region and located approximately in the middle of the c2d area. 
This region contains numerous YSOs. Two notable dark patches show up in the mid-infrared 
Spitzer data, one to the southeast of "Cluster B" and the other to the northeast. Just to 
the south of our JHKg maps is the Herbig Ae star VV Ser surrounded by nebulosity. 



Data Reduction 



We reduced our JHKs data using PyRAF. Each FLAMINGOS field consisted of mul- 
tiple dithers with small offsets around a central position. To reduce our data, we first 
subtracted a median filtered dark image from each dither. Next, we flat-fielded the dithers 
by dividing out a median filtered dome flat. The last, and most important, step was to 
subtract out the sky. This is critical for IR observations where the sky is typically much 
brighter than the astronomical objects of interest. We used a two-pass sky subtraction with 
a median sky computed from the nearest observed dithers in time. We identified the stars 
from the first pass, then masked them out in the second pass to produce an improved median 
sky image. 

The stars in each field were found using the IRAF task daof ind and daophot was used 
for the point spread function (PSF) photometry. We visually inspected the results to remove 
false sources found by daof ind and also to add real sources missed by it. Finally, we used the 
2MASS catalog to correct the coordinates in each field and also to calibrate the photometry. 

The Spitzer data were processed with the c2d data pipeline and are available on the 
Spitzer Science Center's websit^. Details on data reduction may be found in the c2d data 
delivery documentation, also available on the webs i te, and in several previously published 



papers (jPadgett et al.l 120071 : |j0rgensen et al.l l2006l : iRebuU et al.l 120071 : iHarvey et al 
2007bh . 



2006 



3.1. Data Quality 

The mean difference in position between sources in our JHKs catalogs and in the 2MASS 
catalogs is O'.'IS and 95% of our sources have a difference in position of < 0'.'6 (around 2 pixels 
on the FLAMINGOS CCD). To test our photometric error, we computed the difference in 
flux between sources present in our Kitt Peak observations and the 2MASS catalogs. The 
resultant distribution is Gaussian with a = 4 — 5%. 



http : //ssc . spitzer . caltech. edu/legacy/c2dhistory .html 
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The c2d delivery documentation discusses three sources of uncertainty: statistical, sys- 
tematic, and absolute. The first two of these are incorporated into the photometric uncer- 
tainties listed in the c2d catalog. These errors are derived from the repeatability of flux 
measurements using the c2d pipeline. Since we are using the c2d data in this paper, we will 
adopt their uncertainties. The c2d documentation lists a 4.6% systematic error for IRAC 
and 9.2% for MIPS. Combining the systematic and statistical errors, our final photometric 
uncertainties are approximately 5% in IRAC and 10% in MIPS. The absolute uncertainties 
in the flux calibration are 1.5% and 4%, respectively, for the IRAC and MIPSl (24 fim) 
bands. We obtained these values from the Infrared Array Camera (IRAC) Data Handbook, 
Version 3.0 and the Multiband Imaging Photometer for Spitzer (MIPS) Data Handbook, 
version 3.3.0. The absolute uncertainties are added in quadrature with the photometric er- 
rors during source classification, but are not listed in the final catalogs since they are smaller 
than the photometric errors. 

Our 5 and lOcr detection limits are listed in Table [2l For the JHKs bands our lOcj 
limits are 19.5, 18.8, and 17.7 magnitudes, respectively. These hmits are ~ 3.5 magnitudes 
deeper than 2MASS in all three bands. 



3.2. Extinction and Star Classification 



The standard c2d pipeline uses all available wavelengths to compute the extinction 
towards each star: 



log(F,b,(A)/F„,rfeKA)) = log(fc) - 0.4 X Ct{\) x Ay (1) 

where FmodeiW is the stellar photosphere model, k is the scaling factor of the model for 
a particular star, and Cext(A) = A\/Av is the ratio of extinction at wavelength A to visual 
extinction from the dust extinction law. k and Ay are derived from the linear fit of this 
equation by adopting stellar photosphere and dust extinction models. The stellar photo- 
sphere models for /Tg-MIPSl bands are based on the Kurucz-Lejeune models and come from 
the SSC's "Star-Pet" t oofl. For the 2MA SS bands, c2d translated the observed J - H and 



H — K colors of stars (iKoornneefl Il983l ) to fluxes relative to K band and ignored the dif- 
ference between the K and Kg bands. Sources that were a good fit to this equation were 
classified as reddened stars; those that did not fit were compared with other templates to 
classify the sources. 



^http : //ssc ■ spitzer . caltech.edu/tools/ starpet] 
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In this paper, we slightly modified this procedure. First, the extinct ion, Ak,, was com- 



puted for each source using the JHKg bands and the NICER technique (ILombardi &: Alves 



20011 ). We also computed extinctions in A^^ rather than the more traditional Ay so we can 
directly compare our results with those of other authors. The NICER technique relies on 
the assumption of intrinsic values for the J — H and H — Kg colors of stars. For this paper 
we adopted the intrinsic colors J — H = 0.63 ± 0.16 and H — Kg = 0.21 ± 0.14. These values 
were based on the average color for stars from off-cloud fields. Then, using Equation [H we 
identified the stars with the extinction held fixed. 

We made this change because our primary goal in this paper is to compute the extinction 
law at the IRAC and MIPS wavelengths. The method we use for deriving the mid-infrared 
extinction law [Ax/AkJ is dependent on A^^ so we cannot use these wavelengths when 
computing A^^- However, we still need to choose an appropriate extinction law. This will 
be used to compute Ak^ from the JHKs bands and also for the so urce classification. For 



purpo ses of source classification, we considered two dust models, the IWeingartner fc Draine 



( I2OOII ) Rv = 3.1 and Ry = 5.5 models. The Ry = 3.1 model (hereafter WD3.1) is designed 
to reproduce the extinction law of the diffuse interstellar medium; the Ry = 5.5 model 
(hereafter WD5.5) fits the observed law of denser regions. As we discuss in the next section, 
we adopted the WD5.5 dust model for computing the extinctions and source classifications. 
The parameter Ry = Ay/{AB — Ay) was historically used to characterize changes in the 
extinction law. 

The advantage of using only the near-infrared bands to compute A^^ is that the extinc- 
tion law in the near-infrared JHKs bands is very similar for both the WD3.1 and WD5.5 
models. In § l5.3l we will discuss possible biases in our results attributable to our assumption 
about the near-infrared extinction law. 



3.3. High-Reliability Star Catalogs 

We initially created high-reliability star catalogs by selecting only those sources classified 
as stars with both the WD3.1 and WD5.5 dust models. This selection excludes questionable 
stars that can affect our results, but may also miss bona-fide stars that trace a true change 
in the extinction law. To check whether this happens, we plotted in Figure H] those sources 
from Serpens classified as stars using either WD3.1 or WD5.5, but not both. We found that 
the sources classified as stars with WD5.5 are strongly associated with the higher extinction 
regions while the WD3.1 stars are scattered randomly throughout the observed regions. 
Furthermore, up to a third of the detections that are only classified as stars using WD3.1 
have a quality fiag identifying them as confused with a nearby source suggesting that many 
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of these sources have suspect photometry. In the remainder of the catalogs, < 1% of the 
sources are similarly confused. 

For these reasons, we selected all sources identified as stars with the WD5.5 dust model 
and computed line-of-sight extinctions based on this model. We did not want a prominent 
absorption or emission line to bias our results. Therefore, we excluded any sources classified 
as stars only when one of the wavebands was dropped. Furthermore, because accurate 
extinction measurements are essential to this paper, we required sources to have a detection 
> 7a in each of the JHKg bands. Lastly, we 'cleaned' our star catalogs to remove suspected 
faint background galaxies that were mis-identified as stars. This procedure will be discussed 
in detail in §3.51 In our final catalogs we have the following numbers of sources before (after) 
cleaning: 2,403 (2,365) in Ophiuchus, 11,479 (11,280) in Perseus, and 49,712 (49,485) in 
Serpens. 



3.4. Young Stellar Objects and Background Galaxies 



In addition to stars, the c2d classification procedure also identifies Young Stellar Ob- 
ject candidates (YSOc) and background galaxy candidates (Gale). Within the regions we 
observed with deep JHKg and Spitzer observations, we identified the following numbers of 
YSOc (Gale) sources: 25 (6) in Ophiuchus, 40 (25) in Perseus, and 74 (21) in Serpens. Of 
course, many more YSOc and Gale objects are identified in the full cl oud catalogs. Readers 
interested in the YSO/Galc c o ntent of th e full clouds s hould consult Padgett et al.l (2007 ) 



(Ophiu chus), IJorgensen et al.l (l2006l ) and iRebuU et al.l ( 120071 ) (Perseus), and iHarvey et al. 



(l2007al ) (Serpens). All but fourteen of our YSOc sources have the same identification in the 
c2d catalogs. Thirteen of the fourteen are attributable to our modified source classification 
procedure. Our method tends to be conservative when classifying stars, thus ambiguous 
sources may become YSOs. In contrast, the c2d pipeline is more conservative when selecting 
YSOc sources, so ambiguous cases may become stars. The fourteenth source is classified as a 
YSOc in our catalogs, but as a Gale by c2d. This is because it is detected in our deep JHKs 
data but not by the 2MASS survey used by c2d. With our data, the c2d pipeline also clas- 
sifies this source as a YSOc. In the c2d catalogs, this is source SSTc2d J034214. 9-1-314758. 
All of our Gale sources are also identified as such by c2d. 
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3.5. Misidentified Background Galaxies 

Figure [5] is a color-color diagram showing J — H versus H — Kg for our three clouds. 
Those sources brighter than 15th magnitude at Kg are shown in black while fainter sources 
are gray crosses. We expect to see extincted stars following a reddening vector that stretches 
the data points into a line with some scatter, as shown in the figure. We do see this trend with 
the black points in all three clouds and also with the gray crosses in Ophiuchus and Serpens. 
However, this trend is not clear in Perseus where many gray sources with H — Kg ^ 0.7 and 
J — < 1.4 do not appear to follow the reddening vector. 

What is the nature of these sources? One possibility is that these are merely faint stars 
with bad photometry. It is true that these sources are fainter than 15th magnitude at Kg, but 
we are already using our high reliability catalogs to ensure good photometry. Upon further 
examination, we found many of these anomalous sources are only detected in the J-IRAC2 
bands. Another possibility, therefore, is that these sources are not really stars, even though 
they have been classified as such. Perseus is located at Galactic longitude / = 160°, almost 
in the Galactic anti-center direction, while Ophiuchus and Serpens are more towards the 
Galactic center at / = 353° and / = 30°, respectively. With these Galactic longitudes and 
given that these anomalous sources are seen only in Perseus, we suspect these sources are 
likely misidentified background galaxies. Our hypothesis is supported by the recent work of 



Foster et al.l (120081 ). They identified galaxies from near infrared JHKg data and found that 
their galaxies were located in the same region of the J — H vs. H — Kg color-color diagram 
as our suspected galaxies. 

Our goal in this section is to develop a set of criteria that we will use to separate these 
suspected faint background galaxies from our actual stars. It is likely that any procedure 
will remove some bona fide stars and embedded objects in addition to background galaxies. 
We want our criteria to be stringent enough so that the final sample of sources selected 
for removal will be strongly enriched in background galaxies but also contain a minimum 
number of actual stars. We will then use this procedure to create purified star catalogs. 

To determine whether the sources with anomalous colors really are background galaxies, 
we need a group of "known" stars and background galaxies with which we can compare. To 
make these "known" samples, we plotted Kg vs. Kg — [24] for our three cloud regions in 
Figure Ei We only plotted those sources with S/N> 3cr at 24 /im in addition to the earlier 
cut (7(j) we made for uncertainties i n JHKg photometr y. The shaded contours in each panel 



are the SWIRE data of ELAIS Nl (ISurace et al.ll2004l ) as processed by the c2d team. The 



ELAIS Nl region is near the north Galactic pole, and therefore the data should contain 
nothing but stars and background galaxies making them useful for selecting these types of 
objects. Based on Figure Owe selected two groups of sources: those with Kg — [24] < 1 and 
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those sources with Kg > 15 and Kg — [24] > 4. The first group, with no color excess, contains 
stars, while the second population consists of background galaxies. We will use these two 
groups as our "known" samples. 

In Figure [7] we overlaid the "known" stars (white circles) and background galaxies 
(dark gray circles) on our J — H vs. H — Kg color-color diagram. These two classes of 
sources are roughly separated. The stars chiefly lie along the stellar reddening vector, while 
the background galaxies appear below and to the right of this vector, exactly where our 
anomalous sources appear in Perseus. As our first selection criterion, we drew dashed lines 
to separate the stars from the possible galaxies. Note that the separation between the two 
is not perfect, but does establish that sources with J — H > 0.6, H — Kg > 0.6 and to the 
right of the sloped line are more likely to be galaxies rather than stars. The equation for the 
sloped line is: J -H = 1.9 x {H - Kg) - 0.16. 

As noted earlier, approximately three-fourths of the objects in our clouds that lie in the 
region selected by our first criterion are faint enough that they are only detected in the J- 
IRAC2 wavebands. Therefore, in Figure [S] we plotted if — i^s vs. i^s — [3.6] and vs. Kg — [4.5]. 
The gray crosses are those sources selected in all three clouds from Figure [7] by the dashed 
lines. Again, white circles are stars and dark gray circles are background galaxies. In both 
of these plots, the stars and galaxies are clearly separated. We defined a single line to use 
for both plots: H — Kg = 1.32a; where x is either Kg — [3.6] or Kg — [4.5]. 

We combined these two selections to reduce the number of background galaxies in our 
data. Starting from the high-reliability "star" catalogs, we removed all sources satisfying 
the following two criteria: 

1. J -H > 0.6, H -Kg> 0.6, and J -H < 1.9 x {H - Kg) - 0.16 

2. H -Kg< 1.32 X {Kg - [3.6]) and H - Kg < 1.32 x {Kg - [4.5]) 

Our cleaning procedure removed 1.6% of the sources in Ophiuchus, 1.8% in Perseus, 
and 0.5% in Serpens. 



Extinction 



As mentioned in ^ \'d.2l we computed line-o f-sight extinctions, A^^, to each star using 
the NICER technique ( iLombardi &: Alvesll200ll ). These line-of-sight (LOS) Ak^ values give 
us a randomly distributed sampling of the true Ak^ in each cloud. To convert these LOS 
extinctions into a uniform map, we overlaid a grid on top of the data to represent the 
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cells in the final extinction map. At each grid position, we used all of the LOS extinctions 
within a given radius. The extinction value in that cell is then the average of the individual 
extinctions, weighted both by uncertainty and distance from the center of the cell using a 
Gaussian weighting function. 

For each cell we chose an integration radius equal to the full width at half maximum, 
or 2.3548cr. All that remains then is to choose the appropriate resolution for our maps. To 
obtain good statistical accuracy for each cell, we used a resolution that yielded on average 
20+ stars in each cell. For Ophiuchus and Perseus, this corresponds to 90" resolution while 
in Serpens, which is towards the Galactic center, we were able to make a 30" resolution map 
and still have on average 36 stars per cell. 

Foreground stars, which are correctly assigned a small A^^ , bias the average extinction 
for a given cell. We dealt with foreground stars as follows: for each cell with more than 2 
stars, we computed the mean and median Ak^ value. If these two statistics differed by more 
than 25%, then we dropped the source with the lowest Ak^ and recomputed the mean and 
median. However, if the difference between the old and new mean extinction is < 25%, then 
we re-added the dropped source to our catalog. After we identified all the sources to drop, 
we recomputed the extinctions in every cell. Our extinction maps of Ophiuchus, Perseus, 
and Serpens are shown in Figures [9] - [TTJ 

Our Ophiuchus region has the largest dynamic range of A^^ in our clouds, up to 4 
magnitudes. A loop of extinction is located at ~ 16^28'^30^ in addition to the high extinction 
region L1688 on the western side of the map. Two-thirds of the YSOc sources are located 
in this half. 

Our Perseus extinction map contains dust filaments in the upper half of the image 
and little extinction above the A^^ = 0.5 level elsewhere in the map. A cluster of YSOs 
identified as IRAS 03388+3139 is located at RA ~ 3''42"" DEC - 31°47'. Furthermore, the 
high extinction regions appear to be correlated with the darker patches in the color image. 
The dust lane that crosses the nebulous region of IRAS 03382+3145 (~ 3''41'"30" +31°55') 
shows up as a high extinction filament in our extinction map. Almost all the Gale sources 
are in the low-extinction southern half of the map. 

The Serpens extinction map has a higher resolution, 30", than either the Ophiuchus or 
Perseus maps. Two regions show up as holes in our maps because of a lack of background 
stars. The larger hole in the center of the map contains the star-forming region known as 
Cluster B while the extinction hole west of Cluster B shows up as a dark patch even at 
24 yum. 
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5. Dust Properties 

The WD3.1 and WD5.5 models are nearly identical in the near-infrared JHKg bands. 
We will exploit this constancy of the extinction law in the near-infrared by using the ex- 
tinctions computed from the JHKg bands and then extrapolating the extinction law to the 
Spitzer wavebands. 

Starting from our basic equation relating flux and extinction, 

log(F,fc,(A)/Fwd(A)) = log(A;) - 0.4 x ^.^(A) x Ak, (2) 
we can re- arrange this to solve for Cext, defined as Ax/Ak/- 

CeM) = -^[log(A;) - log(F„,,(A)/F^o,ez(A))] (3) 

Since we are only interested in differences in the extinction law, we subtracted CextiKg) 
to eliminate the need for k, the scaling factor. Furthermore, because Cext{Ks) = 1 and 
Fmodei(yKs) = 1 (all stellar models are scaled relative to Ks), the equation simplifies to: 

CextiX) = ^ [\og{FoUKs)/Foi>sW) + logF^odeiiX)] + 1 (4) 
For FmodeiW we use the average stellar model derived in Appendix El 



5.1. in the Clouds 

Our goal in this section is to construct two-dimensional maps of the changes in the dust 
properties across each cloud. We will do this by computing a reduced value along each 
line-of-sight and then applying the same method we used to create Ak^ maps, except this time 
our "extinction" values will be values. We define cis the sum of the difference between 
our computed extinction law and a theoretical model over the Spitzer IRAC wavebands 
(3.6 — 8 /im). The reduced is then simply divided by the number of bands summed 
minus 1, since we have 1 free parameter in the model, namely Ak^'- 



n 



1 ^|^ Qx^^(A)-C-.f-'(A) y 
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^extW is computed from Equation HI C^°'^^^ is the extinction law for a given dust 
model, and ax is the uncertainty in C"^^. We used both the WD3.1 and WD5.5 dust models 
in creating our maps. Note that we only summed over the IRAC bands. As we will see 
later, the extinction law at 24 /xm is almost always significantly higher than either the WD3.1 
or WD5.5 extinction models (§ l5.2.2p . We do not want to bias our results based on those 
data so we excluded this data point when calculating the reduced x^- Finally, we excluded 
any negative values for Cg^j(A) as unphysical. 

We created maps of the values using the same process we used for extinction mapping 
and with the same resolution. Our maps for Ophiuchus, Perseus, and Serpens are shown 
in Figures [12] - [HI We drew our contours starting at = 4 because we observe a definite 
transition between the WD3.1 and WD5.5 dust models at this approximate value (see 
Figure [T5l). Statistically, = 4 would arise by chance about 5% of the time for n = 2. 

In broad terms, high extinction regions are correlated with high Ry = 3.1 x^ values. 
This same correlation is not seen with Ry = 5.5. To make a quantitative comparison between 
X^ and Ak^, we binned our data in A^^ and for each bin determined the average x^ with 
both the WD3.1 and WD5.5 models. Our results are shown in Figure [TH The Ry = 3.1 
X^ curve is shown in black while Ry = 5.5 is shown in gray. At low extinction, Ak^ ^ 1, 
the Ry = 3.1 and 5.5 curves are very similar suggesting it is difficult to distinguish between 
the dust models with our method. Above Ak^ ~ 1, the Ry = 3.1 curve rises sharply while 
Ry = 5.5 does not. 

Our increase in Ry = 3.1 x^ values in denser regions means the WD3.1 model does not 
fit the data in these regions. The WD5.5 model, which includes larger dust grains, produces 
much smaller x^ values, suggesting grain growth is occurring in the denser regions. One 
prominent region within the Perseus Ry = 3.1 x^ map does not fit into this picture. This 
region is located at ~ 3^0.5" +31°30' and circled in red in Fi gure [131 This region has a 
peak x^ value of 39, but since it is non-circular and larger than the integration radius, it 
seems unlikely that a single bad data point is affecting the x^- The puzzling thing about 
this region is that the extinction is fairly unremarkable suggesting that grain growth may 
not be causing the increase in x^- This region will be examined in more detail in § 15.2.31 

5.2. The Mid-Infrared Extinction Law 

We have line-of-sight measurements of Ak^ and the extinction law for every star within 
our clouds. In the previous section, we turned these measurements into maps of the x^ 
deviation from the WD3.1 or WD5.5 extinction laws. Now we will examine the extinction law 
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as a function of wavelength and Ak^ ■ We started by binning our line-of-sight measurements 
into four A^^ ranges from low to high: < A^^ < 0.5, 0.5 < A^^ < 1, 1 < A^^ < 2, and 
Ak^ > 2. For each bin, we excluded any negative C°^l values as unphysical and computed 
the weighted average CextW- If we had used these negative C°^^ values, then our extinction 
law for < Ak, < 0.5 would be Ax/Ak, ~ 0.02 for IRAC2-4 (4.5 - 8 /xm), suggesting 
the dust causes almost no extinction for these wavelengths and far below the predictions of 
dust models. The impact of negative C°^l values on our other extinction bins is negligible. 
Our results for our three clouds may be seen in Figures dSl - [IHl The errorbars in each 
bin represent the minimum uncertainty due to systematic errors in measuring the fluxes, 
propagated through Equation HI 

We plotted three theoretical extinction laws in our figures. The WD Ry = 3.1 and 
WD Ry = 5.5 models are the same WD3.1 and WD5.5 models we used in computing the 
X^. In addition to those two, we plotted a third model, labeled KP v5.0 (Pontoppidian 
et al., in prep). This model is one from a grid of models constructed starting from the 



Weingartner fc Draind (120011 ) parameterization of the grain size distribution. Icy mantles of 
water and other volatiles were then added. The specific model we use from this grid is the 
one with the "best fit" to the c2d mid-infrared extinction law and ice features. Several ice 
absorption features in this model can be seen in the figures, these are due to H2O, CO2, or 
CO. 



5.2.1. The Extinction Law From 3.6 — 8 /im 

In the IRAC bands, the observed extinction laws in all three clouds show a similar 
trend: At low A^^, the extinction law is more consistent with the WD3.1 law while at higher 
extinctions the observed law gradually flattens to become more consistent with WD5.5. In 
Table [3] we list our average extinction law, made from all three of our clouds, along with the 
results from other authors and also our three dust models. 

In the lowest extinction bin, < A^^ < 0.5, the 3.6 — 8 /im extinction law coincides 
with the law predicted by the WD3.1 dust model, except for A^,^/Ak^ = 0.28, which is much 
higher than the WD3.1 value of 0.17. The data from Perseus and Serpens are the cause of 
this; both have an anomalously high A^ s/Ak^ value. If we compute the statistic and 
the associated probability, then these data have a 94% chance of fitting the WD3.1 model 
(;^2 ^ Q 39)^ Q^iy ^ 22% chance of fitting the WD5.5 model (x^ = 4.5). 

The extinction law in the 0.5 < Ak^ < 1 bin appears to be a transition between the 
WD3.1 and WD5.5 models. The WD3.1 is 6.8, only an 8% probability of the data 
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matching the model, while the WD5.5 is 3.8 (29% probability). For 1 < Ak^ < 2, our 
extinction law is very close to the WD5.5 dust model. It is extremely unlikely that the data 
fit the WD3.1 model here (x^ = 80, 2.5 x 10~^^% probability). Compared to this, the WD5.5 
model has = 1-2 with a 76% probability. Lastly, for A^^ > 2, our average extinction law 
is slightly flatter than WD5.5. This is due to the data points from Perseus and Serpens, 
and is primarily because of the 5.8 /xm data point. The presence of water ice may explain 
the 5.8 /xm relat ive extinction. T he KP v5.0 dust model incorporates the 6.02 /xm water ice 
absorption line (IGibb et al.ll2004l ). which is within the IRAC 5.8 fim bandpass. The WD3.1 
model is even more strongly rejected in this bin {x^ = 303, 1.9 x 10^^'^% probability), though 
the WD5.5 model is also a relatively poor fit (x^ = 9.3, 2.5% probability). 

As can be seen in Table [31 our 3.6 — 8 p,m ex t inctio n law for A^, > 1 is consistent 
with the results obtained by other authors. iLuta (Il999l ) measured hydrogen recombina- 
tion lines towards the G alactic center and computed the extinction law from 2.6 — 19 fim. 



Indebetouw et al.l (120051 ) used the IRAC instrument on-board Spitzer and measured the ex- 



tinction l aw along two lines of s ig ht using data fr o m the GLIMPSE Spitzer Legacy Science 
program (IBenjamin et al.l 120031 ). [Flaherty et al.l (120071 ) also used Spitzer to measure the 
extinction law toward five regions. All of these results are in rough agreement; they find a 
fiat extinction law from 3.6 — 8 /xm, consistent with the WD5.5 model. 



5.2.2. The Extinction Law at 24 /im 

At 24 /xm, the extinction law is generally higher than that predicted by any of the three 
dust models. This means more extinction exists for a given column density, or alternatively, 
our observed 24/im fluxes are fainter than models predict. The exact value of ^24/^X3 ranges 
from 0.28 to 1.1 and varies among clouds and extinction bins. The only A24/Ak^ value which 
is not higher than any dust model prediction is for Ak^ > 2 in Ophiuchus, which is only 
marginally higher than the predicted value from WD5.5. However, this value is computed 
from just 2 stars. This general discrepancy is the reason we excluded the 24 /xm data when 
creating our maps in § 15.11 

The value of A24,/Aks appears to decrease as Ak^ increases. This behavior is puzzling 
because if grain growth is occurring in denser regions, as suggested by the 3.6 — 8 /xm 
extinction law, then we would expect regions with a larger column density to have a higher 
value of A24/AK,, not a lower one, since larger grains should have a larger relative extinction. 
However, we observe the opposite trend. 

A likely explanation for this trend is that the average stellar model we are using is 
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incorrect at 24 /im. We can obtain an upper-limit on the value for the average stellar model 
flux at 24 fim. There are 23 stars with 24 fim fluxes with < A^,, < 0.5 in all three clouds. 
If these sources had zero extinction, then the ratio of 24 /xm to Kg band flux for these stars 
would provide an estimate of the true average stellar flux. However, these stars do have 
small amounts of extinction, which will increase the 24 fim to Ks flux ratio, assuming the 
true relative extinction at 24 /xm is less than at Kg. Thus, our ratio of 24 /xm to Kg band 
flux is an upper-limit on the true average stellar flux. For the 23 low-extinction stars in our 
clouds, the flux ratio is 0.014 ± 0.002, the same value listed in Table HI Since the average 
extinction for these stars 0.3 Ak^, the true value will be less than this. As a test, we took 
our combined catalogs for all three clouds and used a smaller average model flux of 0.012 
mJy at 24 fim. In our four extinction bins, our relative extinction becomes (from low to 
high AkJ: 0.65 ± 0.31, 0.56 ± 0.14, 0.52 ± 0.08, and 0.28 ± 0.12. Using a lower value for 
the 24 fim stellar flux almost completely eliminates the trend of ^24/^4^^^ decreasing as A^^ 
increases. The only extinction bin which shows a significant drop in the relative extinction 
is for Aks > 2, but there are only three 24 fim stars in this bin. A larger sample of stars is 
needed to confirm this low value. 

Other authors also found A24:/Ak^ to be higher than predicted by dust models. Their 
values are in rough agreeme nt with ours, as s umin g we use the lower value for the average 



stellar flux proposed above. iFlaherty et al.l (120071 ) were able to measure A24,/Ak^ for two 



of their five regions and found A^/^j A^^ = 0.44 ± 0.02 and 0.52 ± 0.03 for Serpens and 
NGC20 68/71, resp ectively. However, they caution that their sample size is small. Even 



though iLuta (119991 ) did not measure the extinction longward of 19 fim, if their flat extinction 



law was projected out to 24 fim, it would also be ~ 0.5. 

It is clear that the extinction law at 24 fxm is not well understood. With some assump 



tions about the average stellar flux, our findings are consistent with those of iFlaherty et al. 



(120071 ). Both their and our results suggest ^24/^^^ ~ 0.5, or 2 — 3x larger than the value 
predicted by current dust models. New dust models will need to incorporate additional ex- 
tinction at 24 fim to reproduce these observations. Furthermore, both the Flaherty et al. and 
our results are from relatively small samples. A larger study is needed before we can answer 
questions about variations in the 24 fim extinction law with column density or line-of-sight. 



5.2.3. Perseus Map 

In this section we will address the 'anomalous' region in the Perseus Ry = 3.1 map. 
This region is circled in red in Figure [131 53 stars are within the region defined by the 
= 4.0 contour. We plotted the extinction law derived from these 53 stars in Figure [12] 
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and found the extinction law for these data are significantly higher than the WD5.5 model 
with the largest deviation occurring at 5.8yum. The water ice feature at 6.02/im may explain 
the discrepancy, however ices are seen towar ds dense cores arid the average ext inction for 
this region is a relatively modest A^^ = 0.60 (IGibb et al.ll2004l : iKnez et al.ll2005l ). Another 
possibility is that the sources in this region may not really be stars, but a search of the 
NASA Extragalactic Database and Simbad within 2' of 3^^40.5™ +31°30' yielded no matches. 
A final possibility is that the near-infrared extinction law is different in this region compared 
to our assumed WD5.5. We will see how this impacts our mid-infrared extinction law in the 
next section. 



5.3. Changes in the Near-Infrared Extinction Law 

As we mentioned in the introduction, t he near-infra red extinction law can be fit by a 



power law. Ax o: X ^, with f3 = 1.6 — 1.8 (jPraind l2003l . and references therein). For this 



paper we assumed the WD5.5 model, with f3 1.6, is correct in the near-infrared. This 
assumption allows us to directly compare our results with those of other authors, who have 
made the same assumption. Furthermore, our data agree with the /3 = 1.6 reddening vector 
shown in Figure [51 However, if the near-infrared extinction law is actually steeper in some 
regions this will introduce another source of error into our measurements of the mid-infrared 
extinction law. To measure the effect of choosing a different value for (3, we modified the 
WD5.5 law to have f3 = 1.8 in the near-infrared, then re-classified the stars, computed Aks, 
and derived the mid- infrared extinction law. Our values of Ax/Ak^ are decreased by using 
P = 1.8 as opposed to WD5.5. The decrease is ~ 15% for IRACl and IRAC2, ~ 25% for 
IRAC3 and IRAC4, and ~ 15% for MIPSl. The extinction law in Figures [MH] still flattens 
as Axg increases, but does so more gradually. A different value of f3 may partially explain 
the 'anomalous' region in Perseus. Using the (3 = 1.8 law for this region, we found that 
the derived extinction law in the mid-infrared is decreased to the point that it is consistent 
with WD5.5. However, the extinction law is still flatter than what is observed elsewhere in 
Perseus at the same A^^- 

Based on this analysis, we expect our assumption of the WD5.5 extinction law in the 
near-infrared may introduce an error of up to 15% in IRACl, IRAC2, and MIPSl, and 25% 
in IRAC3 and IRAC4. However, we have no way of measuring possible variations in the 
near-infrared extinction law along different lines of sight. Therefore, while it is important to 
keep this in mind as a possible source of error, we have not attempted to incorporate it into 
the errorbars shown in our figures. 
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6. Conclusions 

In this paper we took deep JHKg observations of regions within three molecular clouds: 
Ophiuchus, Perseus, and Serpens. These were combined with Spitzer c2d data to create 
band-merged catalogs in eight wavebands from 1.2 — 24 fim. Our final catalogs contain 2,365 
stars in Ophiuchus, 11,280 in Perseus, and 49,485 in Serpens. From these catalogs we created 
uniform maps of the extinction and x^- We also computed the average extinction law in four 
extinction bins: < Ak, < 0.5, 0.5 < Ak, < 1, 1 < Ak, < 2, and Ak, > 2. 

Our results in the IRAC bands, 3.6 — 8 /im, show that grain growth occurs in the 
dense regions of our c l ouds. The maps show the deviations from two dust models: 



Weingartner fc Draind (120011 ) with Ry = 3.1 and Ry = 5.5 (abbreviated as WD3.1 and 
WD5.5). There is a strong correlation between high extinction and high WD3.1 x^? but not 
with WD5.5 x^- This suggests denser regions are more consistent with the WD5.5 extinction 
model, which has grain sizes up to 10 x larger than those in the WD3.1 dust model. The 
same behavior is seen in the average extinction law. For < A^^ < 0.5, it is most consistent 
with WD3.1 (94% probabihty of the data fitting this model). As Ak^ increases, the extinc- 
tion law gradually flattens and becomes more consistent with WD5.5 for 1 < Ak^ < 2, with 
a 76% probability. For Ak^ > 2, the WD5.5 model is not as good a fit to the data, primarily 
due to the A^^ s/Ak^ relative extinction, which is higher than the model. This may indicate 
the presence of water ice in the densest regions. 

Not all areas in our clouds can be explained using grain growth. Our map of Perseus 
reveals a region that has only a modest average extinction of A^^ = 0.60, yet has a peak 
WD3.1 of 39. The extinction law in this region is much fiatter than WD5.5. Because 
of the low average extinction here, it doesn't seem likely that grain growth is altering the 
extinction law, nor that ices have formed on the dust grains. It is possible the background 
sources in this region are not really stars, or that the near-infrared extinction law is steeper 
in this region. More work is needed to understand the physical cause of the extinction law 
in this region. 

Unlike in the IRAC bands, the 24 fim extinction law is 2 — 4x higher than current dust 
models predict. There is also a negative correlation between Ak^ and the value of A24:/Ak^- 
Both of these results maybe be partially mitigated by considering the small numbers of stars 
detected at 24 /im or slight changes to our average stellar model fiux. However, even with 
the most conservative assumptions, our average value for ^24/^4^^^ is still ~ 2x larger than 
that predicted by du st grain models. This results agrees with the measured values from 



Flaherty et al.l (120071 ). Future dust models will need to reproduce the relatively large value 



of the 24 /im extinction. 
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A. Average Stellar Models 

To compute th e ex pected stellar distrib ution for each cloud, we used the galaxy model 



from I JarrettI (jl992l ) and lJarrett et al.l (jl994l ). Given the Galactic coordinates, distance, and 
average extinction for a cloud, this model produces source counts broken down by spectral 
type. These source counts provide a rudimentary weighting function that could be used in 
computing the average stellar model. However, a weighting function from the raw source 
counts would ignore the reality that not all spectral types are equally observable at each 
wavelength. Our next step was to correct for the detection limits of our actual observations. 
The galaxy model lists the extincted Kg magnitude for each star and its extinction. Ay. 
We used this information, the Weingartner & Draine Ry = 3.1 extinction law, and the 
stellar templates from "star-pet"|f| to compute the magnitude of each star at every observed 
waveband from J — 24 /xm. We used WD3.1 since this was the law used by the galaxy model 
when computing the extinctions. We arbitrarily discarded any detection falling below the Sex 
cutoff for the given waveband (Table [2]). Using the normalized, final, corrected source counts 
at each waveband as a weighting function, we then computed the average stellar template 
for each region and an overall average stellar model. Our results are shown in Table HI We 
repeated the above procedure using the Weingartner & Draine Ry = 5.5 extinction law. The 
average fluxes using this law never vary by more than 0.001 mJy from those listed in Table 

SI 



http : //ssc ■ spitzer . caltech. edu/tools/starpet| 
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Table 1. Basic Properties of the Clouds 



Cloud 


I 

(deg.) 


b 

(deg.) 


Dist. 

(pc) 


Area Mapped in JHKg 
(deg^) 


Ophiuchus 


353 


16 


125 ±25^^ 


0.15 


Perseus 


160 


-19 


250 ± 50"^ 


1.0 


Serpens 


30 


5 


260 ± 10" 


0.33 



de Geus et all fll989f ) 



fenoch et all (120061 ) 



^Straizvs et al. (1996) 



Table 2. Magnitude Limits 



Band 


lOo" Limits 


5(7 Limits 


J 


19.5 


20.3 


H 


18.8 


19.5 


Ks 


17.7 


18.4 


IRACl 


16.6 


17.3 


IRAC2 


15.9 


16.6 


1RAC3 


13.6 


14.3 


IRAC4 


12.8 


13.5 


MIPSl 


7.7 


8.4 



Table 3. Average Relative Extinction A\/Ak^ 



Source 




3.6 j2m 


4.5 /im 


5.8 /im 


8 /im 


24 /im 


Our data, < Ak^ < 0.5 




0.41 ±0.19 


0.26 ±0.18 


0.28 ±0.18 


0.21 ±0.17 


1.08 ±0.32 


Our data, 0.5 < Ak, < 1 




0.49 ±0.10 


0.35 ±0.10 


0.35 ±0.10 


0.35 ±0.10 


0.75 ±0.14 


Our data, 1 < Ak^ < 2 




0.60 ±0.05 


0.46 ±0.05 


0.44 ±0.05 


0.43 ±0.05 


0.61 ±0.08 


Our data, Ak^ > 2 




0.64 ±0.03 


0.53 ±0.03 


0.46 ±0.03 


0.45 ±0.03 


0.34 ±0.13 


Other Authors 


Flahertv et al. f2007) 




0.632 ±0.004 


0.54 ±0.01 


0.50 ±0.02 


0.50 ±0.01 


0.46 ± 0.04 


Indebetouw et al. f2005) 




0.56 ±0.06 


0.43 ±0.08 


0.43 ±0.10 


0.43 ±0.10 




Lutz (19mY 




0.53 ±0.03 


0.50 ±0.08 


0.49 ±0.06 


0.42 ± 0.06 




Dust Models^ 


Weineartner & Draine (2001). 


= 3.1 


0.40 


0.25 


0.17 


0.22 


0.17 


Weineartner & Draine f2001'). 


= 5.5 


0.60 


0.49 


0.40 


0.41 


0.24 


KP, v5.0 




0.48 


0.38 


0.34 


0.38 


0.38 



'^The extinction for the closest ISO wavelength to each Spitzer band is listed: 3.7, 4.4, 5.9, and 7.5 /im 
'^Extinctions are computed at the central wavelength in each band 



Table 4. Average Stellar Model for Each Cloud 



Region 


J 


H 




IRACl 


IRAC2 


IRAC3 


IRAC4 


MIPSl 




(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


Ophiuchus 


1.374 ±0.291 


1.335 ±0.095 


1. 


0.448 ±0.018 


0.267 ±0.011 


0.178 ±0.008 


0.103 ±0.006 


0.014 ±0.001 


Perseus 


1.238 ±0.276 


1.284 ±0.099 


1. 


0.457 ±0.023 


0.269 ±0.011 


0.180 ±0.013 


0.101 ±0.008 


0.014 ±0.000 


Serpens 


1.560 ±0.232 


1.393 ±0.061 


1. 


0.443 ±0.013 


0.261 ±0.016 


0.177 ±0.005 


0.106 ±0.005 


0.014 ±0.000 


Average 


1.412 ±0.152 


1.357 ±0.046 


1. 


0.447 ±0.010 


0.267 ± 0.007 


0.178 ± 0.004 


0.104 ±0.003 


0.014 ±0.001 



I 

W 

I 



-24- 




Fig. 1. — Color image of the Ophiuchus cloud made from the 3.6 /im (blue), 8.0 /im (green), 
and 24 fim (red) channels on Spitzer. The green outline denotes the region we observed with 
our JHK^ observations. 
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Fig. 2. — Color image of the Perseus cloud made from the 3.6 fj,m (blue), 8.0 /im (green), 
and 24 fan (red) channels on Spitzer. The green outline denotes the region we observed with 
our JHKs observations. 
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Fig. 3. — Color image of the Serpens cloud made from the 3.6 /xm (blue), 8.0 fim (green), 
and 24 fim (red) channels on Spitzer. The green outline denotes the region we observed with 
our JHKs observations. 
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Fig. 4. — Extinction map of Serpens (See §|1]) with sources identified as stars using only the 
WD3.1 (white) or WD5.5 (dark gray) dust models. 
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h + t<- + .1111-1 



012 012 012 

Fig. 5. — Near- infrared color-color diagrams for Ophiuchus, Perseus, and Serpens. Gray 
crosses have magnitude > 15 while the black points are brighter. The extinction vector 
shown is for the WD5.5 model. 
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Fig. 6. — Ks vs. Ks — [24] plot for our three clouds. We used four different symbols 
to correspond to different source classifications: stars are circles, background galaxies are 
squares, 'YSO candidates' are triangles, and plus signs denote all other classifications. Th e 
shaded contours are from the c2d processed catalog of part of ELAIS Nl (ISurace et al.ll2004j ). 
The dashed lines denote the cutoffs we used for selecting "known" stars and background 
galaxies. 
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Fig. 7. — Near-infrared color-color diagrams for Ophiuchus, Perseus, and Serpens. The white 
and dark gray circles are the 'known' stars and background galaxies, respectively. The dashed 
lines select those sources with J-H > 0.6, H-Ks > 0.6, and J-H < 1.9 x (H-Ks) -0.16. 
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Fig. 8. — H — Kg vs. Ks — [3-6] (left) and vs. Kg — [4.5] (right). The gray crosses are the 
sources selected by the dashed lines in Figure [7] for Ophiuchus, Perseus, and Serpens. The 
white and dark gray circles are the 'known' stars and background galaxies, respectively. The 
dashed line in each panel has the equation H — Kg = 1.32a; where x is either Kg — [3.6] or 
Kg - [4.5]. 
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Fig. 9. — niap of our observed region within the Ophiuchus cloud. The map has a 
resolution of 90". The contours start at A^^ = 0.5 in steps of 0.25 (5cr). We plotted all the 
YSOc and Gale objects as triangles and squares, respectively. The maximum Ak^ value in 
the map is 4.0 magnitudes. 
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Fig. 10. — Ak, map of our observed region within the Perseus cloud. The map has a 
resolution of 90". The contours start at A^^ — 0.5 in steps of 0.25 {5a). We plotted all the 
YSOc and Gale objects as triangles and squares, respectively. The maximum A^^ value in 
the map is 2.6 magnitudes. 
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Fig. 11. — map of our observed region within the Serpens cloud. The "holes" in the 
map correspond to regions with no data eitherbecause they are too dense or because the only 
sources present were not classified as stars. The map has a resolution of 30". The contours 
start at Ak^ = 0.5 in steps of 0.25 (5cr). We plotted all the YSOc and Gale objects for 
Serpens as triangles and squares, respectively. The maximum Ak^ value in the map is 3.4 
magnitudes. 
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Fig. 12. — Map of the Ak^ and values in Ophiuchus. The top panel shows the extinction 
map with contours starting at 0.5 magnitudes in steps of 0.25 mag {5a). The middle and 
bottom panels show the maps for the same region made assuming either the WD3.1 
(middle) or WD5.5 (bottom) extinction laws. Contours start at = 4 in steps of 4. 
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Fig. 13. — Map of the A^^ and values in Perseus. The left panel shows the extinction map 
with contours starting at 0.5 magnitudes in steps of 0.25 mag. {5a). The middle and right 
panels show the maps for the same region made assuming either the WD3.1 (middle) or 
WD5.5 (right) extinction laws. Contours start at = 4 in steps of 4. The region circled in 
red is discussed in §§ 15.11 and 15.2.31 
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Fig. 14. — Map of the A^s and values in Serpens. The top panel shows the extinction map 
with contours starting at 0.5 magnitudes in steps of 0.25 mag. (5cr). The middle and bottom 
panels are the maps for the same region made assuming either the WD3.1 (middle) or 
WD5.5 (bottom) extinction laws. Contours start at = 4 in steps of 4. 
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Fig. 15. — versus Ak^ for the cfouds we observed. The WD3.1 values are black, while 
the WD5.5 values are light gray. The error bars shown are the standard deviations of the 
mean value in each bin. 
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Fig. 16. — The extinction law in Ophiuchus in four different Ak^ ranges: 0-0.5, 0.5-1, 1-2, 
> 2. The data points are the weighted average observed extinction law computed from all 
sources within the specified A^^ range where the errorbars are the minimum uncertainty 
due to systemat ic errors in measuring t he flu xes. We plotted three different dust models 
for comparison: IWeingartner fc Draind (120011 ) Ry = 3.1 and 5.5 ( dotte d, gray lines), and 
KP v5.0 (Pontoppidian et al., in prep) in black. Using iGibb et al.l (120041 ) as a reference, we 
identified the ice features in KP v5.0. 
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17. — The extinction law in Perseus in four different Ak^ ranges: 0-0.5, 0.5-1, 1-2, > 2 
Figure [16] for a description of tlie data points, curves, and errorbars. 
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18. — Tfie extinction law in Serpens in four different Ak^ ranges: 0-0.5, 0.5-1, 1-2, > 2 
Figure [H] for a description of the data points, curves, and errorbars. 
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Fig. 19. — The average extinction law obtained for the region in Perseus with anomalously 
high in the Ry = 3.1 map (see §§ 15.11 and I5.2.3p . The coordinates of this region are 
(~ 3^^40.5™ +31°30'). See Figure[T6]for a description of the data points, curves, and errorbars. 



